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Method for synchronization in a spread spectrum receiver 
BACKGROUND OF THE INVENTION 

5 

1. Field of the Invention 

The invention relates to synchronization methods for spread spectrum receivers. 
Especially, the invention is related to such a method as specified in the preamble of 
10 claim 1. 

2. Description of Related Art 

Wideband Code Division Multiple Access (WCDMA or CDMA) technology is a 
15 strong candidate for future global wireless mobile communications. WCDMA [1] 
has been selected as an air interface solution e.g. in UMTS (Universal Mobile 
Telecommunication System) standard, which will provide a multitude of services, 
especially multimedia and high bit rate packet data. 

20 CDMA is based on spread sprectum technique, where the idea is to spread the 
narrowband information signal into a common wide frequency band before 
transmission. In Direct Sequence (DS) CDMA, the spreading is performed by 
wideband noise-like signal, which simultaneously identifies each user in the system. 
This pseudo-noise is also called the code or the chip sequence of a particular user. 

25 

The final objective in the reception of a DS-CDMA system is to estimate the 
symbols which carry the data, but a prerequisite task is to get the local code 
generator synchronized to that of received signal. This means estimation of the 
propagation delay, which gives the required knowledge to the receiver about the 
30 phase of the spreading code. In addition, the strengths of possible multipaths, and 
carrier phase must be estimated. Code timing estimation tends to be the most 
challenging task, and the other parameters can be estimated given a reliable delay 
estimate [2,3]. 

35 Conventional CDMA systems rely on single-user techniques, such as matched filter 
[41, in both delay estimation and detection. Although simple, they are inadequate if 
the code orthogonality conditions are perturbed. This happens even in synchronous 
system with orthogonal codes due to the existence of multipaths with different 
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delays. Moreover, if the desired signal is much weaker than the interfering signals, 
which is commonly known as a near-far problem, single-user techniques can totally 
collapse. In the uplink (e.g. mobile to base) communication, near-far problem can be 
mitigated by power control, but if used in downlink (e.g. base to mobile) 
5 communications, it actually causes near-far effect for the downlink receiver. This is 
the case e.g. in the WCDMA concept, where equal performance is offered to each 
mobile user, regardless of their locations. The work of Verdu [5] meant, however, 
that it is possible to eliminate the effects of multiple access interference and near-far 
problem. This optimum multiuser detector, although computationally demanding 
10 and requires all the system parameters to be known, initiated the development of 
lower complexity near-far resistant techniques for both detection [6,7] and 
delay/channel estimation [8,9]. 

Code timing acquisition refers to coarse delay estimation, where the maximal 

15 acceptable error is half a chip duration. In [8,9] a subspace approach was 
considered, which was based on multiple signal classification (MUSIC) [12]. 
However, their performance is found to be inadequate in highly loaded systems [10], 
when the signal might not have subspace structure anymore. Since acquisition takes 
place before actual data transmission, training symbols or preamble can be used. 

20 This was utilized in [10] and [11]. Both maximum likelihood -based algorithms 
modeled the training symbols as desired signal and all the interference as colored 
non-Gaussian noise that is uncorrelated with the desired signal. The algorithm in 
[10] was found to be near-far resistant and tolerate high system loading. In addition, 
the accuracy of the delay estimate was achieved by solving a second-order 

25 polynomial for each chip interval. A major limitation, however, is that the method 
can not be extended to fading channels. This is mainly because of the need of long 
training period, during which the channel may have a (sample) mean too close to 
zero. This disadvantage was removed by receiver diversity in [13], which allowed 
shorter training. The algorithm in [1 1] also has limitations in extension to the fading 

30 channels, and moreover looses its near-far resistance in highly loaded systems. The 
computational complexity of all those algorithms is also still quite high from the 
downlink signal processing point of view. Most promising performance with low 
computation was found in [14], where the effects of interference were suppressed by 
correlating the matched filter output with a delayed version of it. The motivation for 

35 this differential correlation is roughly speaking as follows: when only +rs are sent 
to desired user, its contribution in the received signal at the symbol level is only the 
fading process, whose time-correlatedness can be exploited. Basically the same idea 
was exploited also in [15], with the exception that the differential correlation was 
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performed prior to matched filtering. This was due to the special shift-and-add - 
property of m-sequences, which makes the proposed algorithm sensitive to the 
choice of spreading code, 

5 SUMMARY OF THE INVENTION 

An object of the invention is to realize a synchronization method, which alleviates 
the problems associated with the prior art. A further object of the invention is to 
realize a synchronization method, which is more interference tolerant than the prior 
10 art methods. 

The objects are reached by exploiting differential correlation. 

The method according to the invention is characterized by that, which is specified in 
15 the characterizing part of the independent method claim. The dependent claims 
describe further advantageous embodiments of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 Various embodiments of the invention will be described in detail below, by way of 
example only, with reference to the accompanying drawings, of which 

Figure 1 illustrates a channel model used in the description of an advantageous 
embodiment of the invention. 

Figure 2 illustrates a part of a block diagram of a spread spectrum receiver based 
on non-coherent matched filtering according to prior art and a part of a block 
diagram of a spread spectrum receiver based on non-coherent differentially 
correlated matched filtering according to the invention, 

Figure 3 illustrates an example of a typical fading process, 

Figure 4 illustrates the results of a comparison simulation of probability of 
acquisition with inventive and prior art synchronization methods, 
35 

Figure 5 illustrates the results of a further comparison simulation of probability of 
acquisition with inventive and prior art synchronization methods, 



25 



30 
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Figure 6 illustrates the results of a comparison simulation of mean acquisition 
times with inventive and prior art synchronization methods, 

Figure? illustrates the results of a further comparison simulation of mean 
5 acquisition times with inventive and prior art synchronization methods, 

Figure 8 illustrates the results of a comparison simulation of number of acquired 
paths with inventive and prior art synchronization methods, 

10 Figure 9 illustrates the results of a comparison simulation of probability of 
acquisition as a function of number of users with inventive and prior art 
synchronization methods, 

Figure 10 illustrates the results of a comparison simulation of mean acquisition 
15 time as a function of number of users with inventive and prior art synchronization 
methods. 

Figure 1 1 illustrates the results of a comparison simulation of probability of 
acquisition as a function of number of paths with inventive and prior art 
20 synchronization methods, 

Figure 12 illustrates the results of a comparison simulation of RMSE as a function 
of number of users with inventive and prior art synchronization methods, 

25 Figure 13 illustrates the results of a comparison simulation of RMSE as a function 
of preamble length with inventive and prior art synchronization methods, 

Figure 14 illustrates a block diagram of a part of a receiver according to an 
advantageous embodiment of the invention, and 

30 

Figure 15 illustrates a block diagram of a part of a receiver according to a further 
advantageous embodiment of the invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 



Accordiiig to a further aspect of the invention, simple and efficient algorithms for 
the code timing acquisition in the Direct-Sequence Code-Division Multiple Access 
(DS-CDMA) communication system. The essential assumption is that a preamble is 
available for the desired user. Then the correlation matrix R(r) of the sampled data, 
where r is suitably chosen time lag, contains the timing information only of desired 
user, while the contributions of uncorrelated interferers and noise are suppressed 
out. Coarse delay estimates are then obtained by matched filter (MP) or multiple 
signal classification (MUSIC)-type approaches. In the latter case only L eigenvec- 
tors are computed, where L is the number of resolvable paths. For cases in which 
only one path exists, an additional procedure is proposed to both approaches, by 
which the estimation accuracy is greatly improved with neglible increase in compu- 
tation. More precisely, the chip timing offset due to chip-asyncronous sampling can 
be determined by solving a system of two second-order polynomials for each chip in- 
terval. Therefore, only at most 2C hypotheses are needed, where C is the processing 
gain. All the proposed methods are computationally quite simple, containing mainly 
MF-operations, or at most computation of only few eigenvectors. Numerical experi- 
ments speaks for the possibility of achieving significant performance gains compared 
to conventional acquisition, making them attractive options to be attached for the 
next generation mobile receivers. 

The present aspect of the invention brings several improvements over [14]. Accord- 
ing to the present aspect of the invention, we assume periodic (or "short" ) spreading 
codes, and analyze the properties of the correlation matrix of the received and sam- 
pled data with non-zero time lags. In short, we call it a diflferential correlation 
matrix. In addition, two extentions are made. First, the assumption of periodic 
codes makes it possible to consider also eigenvalue- based approaches. We do this by 
applying MUSIC (12]. which has also earlier been applied to CDMA synchroniza- 
tion in [8, 9]. In those papers, however, MUSIC was applied to the data correlation 
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matrix with zero time lag, and problems arose in high system loads due to the dis- 
appearance of signal subspace structure. According to the present aspect of the 
invention, it is reasonable to apply MUSIC even in highly loaded systems, because 
differential correlations effectively filters interference prior to actual delay estima- 
tion. Therefore, the signal subspace dimension is significantly reduced while still 
preserving the desired information, making the circumstances favorable for MUSIC. 
Second, we assume arbitrary delays, and propose a simple way to accurately esti- 
mate the chip timing offset in single-path case. This is an important task since the 
estimation errors in chip timings may cause serious performance losses in symbol 
estimation [16]. The invention enables determining the fractional part of the delay 
by solving a system of two second-order polynomials for each chip interval. 

1. PROBLEM FORMULATION 

1.1 NOTATIONS 

Throughout the paper, lower and upper case boldface letters denote vectors and 
matrices, respectively. In addition, we denote 

(•)* complex conjugate 

(•)^,( )^ transpose, and Hermitian transpose 

E{'} expected value 

sign{ * } sign-function 

x„ nth vector 

II • II Frobenius norm 

0, 1 matrices with all entries 0 and 1, respectively 

1.2 SIGNAL MODEL 

The signal model studied in this paper is a baseband downlink model with fading 
multipath channel, and additive white gaussian noise. The channel model is depicted 
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in Figure 1, where the signal sent by the base station is x{t) = Em=i ^k=i ^kmSkit - 
mT), containing the information of M symbols of K users. Here bkm is kth user's 
mth symbol, Ski ) is Jfcth user's chip sequence. Sk{t) e {-1, +1}, * G [0, T), Sk{t) = 0, 
t ^ [0,T), where T is the symbol duration. We assume that the channel is fixed 
during one symbol, i.e. a<(t) = aim,t € [mT, (m + 1)T) in the Figure 1. aim is 
called an attenuation factor of the Zth path, which a complex number and may vary 
symbol by symbol. The received signal hence have the form 

K*) = E E &ibm E aimSkit -mT- xiTc) + n(t), (1) 

where L is the number of resolvable paths, Tc is the chip duration, and x«^c is the 
delay of the Zth path, where Xi = di + Si with di integer and 6i G [0, 1). The delay of 
each path is assumed to remain roughly constant in the observation interval. n{t) 
denotes noise, and the chip sequence length (i.e. processing gain) is C = ^. From 
now on, we assume Tc = 1 for simplicity. 

First, the received data is sampled by chip-matched filtering, and the equispaced 
samples are collected into C-vectors r„i, 

= [r[mC], r[mC + 1], • • • i r(mC + C - 1]]^. (2) 

They have the well known form (e.g. [8]) 

K L 

Tm = E El^^"»-l^^."*-l£w + aimbkmCkl] + ^my (3) 

where Um denotes noise vector, and 

cjti = ^lixi) = (1 - Si)gj^i{di) + Stg^^idi + 1) (4) 
Cki = Ckiixi) = (1 - + SigM + 1). (5) 

Here the "early" and "late" parts of the code vectors are 

K,Mt)=[sk[C-di + l] ... s,[C] 0 ... o]^ (6) 
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g«(rfO=[o .-.0 Sk[l] Sk[C-di]]^^ (7) 

Notice that Eqs. (4)-{5) and (6)-(7) are due to chip and symbol asyncronous sam- 
pling, respectively. To get more compact representation of the data, define a C x 2KL 
dimensional code matrix G, 

G =^ [£ii,Cii,...,Ci£^,Cix„. ..,Cj<'£,,Eicx,], (8) 

and a 2- vector 

Stacking all vectors Zkim into a 2iirL- vector a„i, 

am = [Znm. • • • . ZUm> • • • i Z^LmJ^i (10) 

the sample vector (3) can be rewritten as 

rm = Ga^ + n^. (11) 

Here G depends on the codes and the delays, and a^n depends on the channel 
coefficients and the symbols. 

2. SYNCHRONIZATION ALGORITHMS 

Our goal is to compute the correlation matrix R(r) of the sampled data (11), where 
the time lag equals r ^ 0 symbols. In short, we call R(r) a diflFerential- instead of 
autocorrelation matrix to stress the nonzero time lag. The proposed methods are 
based on the observation that R(t) contains timing information only of the desired 
user, which is first shown. The following realistic assumptions are made: 



#1. A constant preamble = l,m = 1, . . . , M ("all ones") is available for the 

desired user, which is the user A; = 1 from now on. 
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#2. Symbols of interfering users are uncorrelated random binary variables of zero 
mean. Then E{bkmbk:m+r} = 0, and E{bkm} = 0 for 51^ 1. 

#3. Fading processes of each path are uncorrelated and zero mean. Then E{aimCiiin} = 
0 for / # and E{oi„} = 0. 

#4. Fading process of each path is stationary. Therefore, the correlation coefficient 
equals 

"'(^) = inS^ = r = m - n. and = E{M. 

#5. Noise is a random zero mean variable, and independent with the data and 
fading. Then £;{n„.n^[+J = 0 = E{ 

^m^^+r}) where iim and are those of 

Eq. (11). 

Notice that assumptions (#3) and (#4), without zero-mean assumption, are equiva- 
lent to the common uncorrelated scatterer (US) model [4], and wide-sense stationary 
(WSS) model [18], respectively. It should be noted, however, that assumption (#3) 
is not necessary for the purposes of the paper, but it is made to achieve a simple 
and elegant representation of R(r). The case where (#3) do not apply is discussed 
in Remark 1. 

According to Eq. (11), and assumption (#5), the differential correlation matrix now 
is 

Kir)''^ E{TmV^^,} = £{(Ga^ + n„,)(Ga^+, + n^+,)'^} 

= GA(r)G^, (12) 

where the 2KL x 2KL dimensional matrix A(r) is defined as 

A(r) £;{a,„a^^ J. (13) 
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According to Eq. (10), the matrix A(r) contains {KVf submatrices of dimension 
2x2. Denote these by A(fcj),(jk»j')(r). i.e. 



(14) 



From Eq. (9) we directly have 



^/,m-l^ir',m+T-l^^,m-l6A;',m+T-l 



}• 

(15) 



Since fading and symbol processes are independent, we have moreover 



A(fcj),(fc'j')(r) = 



E{ai,m-l a?' ,m+r- 1 } -E" {fc*.m- 1 h' ,m+T- 1 } 

^{aimaj.,„,+^_i ,m+T-i } 



■E^{ai,m- 1 ai"' .m+T } ^fe.m- 1 6ik' .m+T } 

(16) 



We next show that almost all matrices ^{ki)Xkn*){^) are equal to 0, provided the 
time difference is equal to r = 2 symbols or more. The reason why r = 1 is not a 
good choice for the time lag is discussed in Remark 2. 

Case 1 : 1 = = A;'. Now all the symbols are equal to 1, and A(i/),(if/)(r) depends 
only on the channel coefficients. With I V we have A(H)^(i//)(r) = 0 due to 
assumption (#3). On the other hand, with / = we have due to assumption (#4) 



A(u).(io(^) = /^? 



(17) 



a/(r) ai(r + l) 
a/(r - 1) OLiir) 

Case 2 : 1 = A; 9^ A:'. Now ^{\\)ivv){j) depends on the channel coefficients and the 
symbols of user k\ Since these symbols are zero mean, A(H)^(jt'/')(7') = 0 for all 
Case 3 : 1 7^ /c, A;'. In this case, we always have ^{ki)j<jkfv)ij) = 0 due to assumption 
(#2). 



Hence, the matrix A(r) contains only L non-zero submatrices of 2 x 2. and they 
correspond to situations where A: = A:' = 1 and I = Due to the definition of 
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Eq. (10), these submatrices are situated in the main diagonal of A(r). Bq. (12) 
therefore reduces to 

A(n),(ii)(r) 0 



R(r) = G 



(r) 



(18) 



from which we see that R(r) depends only on 2L first rows of the codematrix 
G. Since these rows correspond to the desired user A: = 1 only, see Eq. (8), the 
contribution of interference is suppressed from R(r). To express in detail the form of 
R(r), we make a mild assumption that a/(r - 1) ?^ ai(r) ^ai{r + l), and denote it 
Q:i(r). This means that the fading process has almost the same degree of correlation 
whether the time difference is r — l,r or r + 1 symbols. Then, Eq.(17) reduces to 
A(i(),(u)(t) = fijai{T)ly and consequently (see details in the Appendix A) 

filaxir)! 0 



R(r) = G 



(19) 



where 



Cii == £if + Cii. 



(20) 



From the structure of R(r) we see explicitly that it contains the code and delay 
information of the desired user A: = 1 only. 

Remark 1. The reason for assumption (#3) was to achieve a compact represen- 
tation Eq. (19). However, to ensure that the differential correlation matrix R(r), 
with T > 2 symbols, contains only the desired user's information, assumption (#3) 
is not necessary. This is because the contribution of interfering users are got only 
via "Case 2" and "Case 3" above. However, in those cases the zero submatrices were 



wo 00/49720 



PCT/FI00/00n9 



12 

always due to the uncorrelatedness and zero-mean properties of the interfering user's 
s\-nibols. Therefore, R(t) is free of contributions of interfering users even without 
assumption (#3), although the form of R(r) is little bit different. This is because 
now in "Case 1", both subcases "/ = and "/ # I'" cause non-zero submatrices. 
Assumption (#3) is valid e.g. for Playleigh-fading channels [4], but do not apply e.g. 
for Rician fading (4) and fixed multipath channels. In the former case, the paths 
are uncorrelated but non-zero mean, whereas in the latter case, the paths are both 
correlated and non-zero mean. 

Remark 2. It is important to notice, that the situation r = 1 is not as desired 
as the case t > 2 symbols. The reason is as follows: due to symbol asyncronous 
sampling, the sample vectors im usually contain information about two successive 
s>-mbols, see Eq. (3). Therefore, with the time lag r = 1 symbol, two successive 
sample vector have some degree of correlation also with respect to the interfering 
users. To see this in detai, let r = 1, and consider 



A(ikj),(fc'r)(l) = 



E{aimaJ,^}E{bkmh'm} 



E{ai,m-l0.l>,m+\}^{^k,m-lbk',m+l} 



(21) 



If now 1 7^ A; = Jk' with / = /', which belongs to the "Case 3" above, we would have 
a non-zero submatrix equal to 



A(w),(w)(l) = 



(22) 



0 0 
E{\aim?} 0 

If R(l) is now rewritten in the form of Eq. (19), R(l) would clearly have a contri- 
bution which corresponds to the interfering user # 1. For this reason, we consider 
only cases where t > 2 symbols. 

Remark 3. It is worth mentioning that the differential correlation matrix R(r) 
with r > 2 contains two kinds of information about the channel: the average powers 
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of the paths, and the rapidity of their fading processes. The latter is related to the 
correlation coefficient ai(r), which is the closer to 0 (resp. 1) the faster (resp. slower) 
is the facing process. Therefore, provided equienergy paths, the slowly varying paths 
are the best represented in R(r). This is a nice property, because such paths are 
more valuable to be acquired due to their easier tracking. 

2.1 MF-TYPE APPROACH 

In practice R(r) is estimated by R(r) = Em=i ^mr^+r- The simplest way to 
estimate the delay is try to match the known code to the data as well as possible, 
i.e. to find the solution for 

X-argmax ^^^^^^^^^^ (23) 

Here Ci (x) is a replica of the user's A; = 1 code, corresponding to the fractional 
testdelay x = + <J- Formally, 

ci(x) = (l-<J)gi(d) + <Jgi(rf + l) (24) 
= (1-5)' 



Si[C-d+l] ••• 5i[C] 5i(l] ••• si[C-d] 



+ 5 



si[C-d] ••• si[C] si[l] ••• si[C-d-l] 



(25) 



The delay estimator Eq. (23), which we label DC-MF (Differential Correlations 
based MF), was the main contribution of our earlier work [19]. Notice that as "all 
ones" preamble is used, the conventional non-coherent MF is exactly DC-MF with 
r = 0. It is also worth to notice that, since Eq. (23) equals 

DC-MF includes, indeed, only MF and delay operations. Estimator of form Eq. 
(26) would be the practical implementation of DC-MF, and is depicted in Figure 2. 
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2.2 MUSIC-TYPE APPROACH 

In practice, the estimate R(r) will always contain some contributions of interfering 
users and noise due to the finite length of the preamble. But if R(r) would have 
exactly the form of Eq. (12), or equally Eq. (19), it could be rewritten by eigenvalue 
decomposition as 

R(r) = U(r)AU(r)^. (27) 

Here A is a diagonal matrix, whose diagonal elements are L non-zero eigenvalues of 
R(r), and U(r) is a C7 x L dimensional matrix, containing the related eigenvectors 
on its columns. But then MUSIC (Multiple Signal Classification) [12] could be 
used to estimate the delays. Due to finite number of samples, R(r) will be in 
practice asymmetric, and it must be first symmetrized to achieve real eigenvalues and 
orthonormal eigenvectors. Therefore, we define symmetric Rsym('r) = R(r)+R(r)^, 
for which MUSIC is applied. More precisely, the delay estimation is performed as 

iic.(x)r ■ 

where U(r) consists of L principal eigenvectors of Rsyni(r). We label this estimator 
as DC-MUSIC (Differential Correlations based MUSIC). 

MUSIC has been applied to CDMA chip timing acquisition also earlier in [8, 9], but 
it should be stressed that DC-MUSIC is totally different from those. This is because 
in those papers MUSIC was applied to the data correlation matrix R(0). Formally, 
as the preamble is "all ones"-type, the delays are estimated by MUSIC according to 

X-argmax . (29) 

where now U(0) consists of N principal eigenvectors of R(0). N is called the sig- 
nal subspace dimension, which is equal to the rank of G, see Eq. (8), or C at the 
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maximum. Here the columns of U(0) are said to span the so called signal subspace, 
which is orthogonal to the noise subspace spanned by the rest eigenvectors of R(0). 
Dimensions of the spaces are min{7V, C} and max{C - N, 0}, respectively. Consid- 
ering the data model (3), N can be as high as = 2KL due to multipath fading, 
and even in non-fading case N is still as large as AT = 2K. In highly loaded system, 
when N > C, R(0) do not obey subspace structure anymore, and this causes the 
failure of MUSIC. In DC-MUSIC, on the other hand, the first thing to do is to filter 
interference and noise out as much as possible by differential correlations, i.e. by 
estimating R(r) instead of R(0). After this MUSIC is applied by estimating the L 
principal eigenvectors of symmetrized RsymC''')* 

TIMING-OFFSET ESTIMATION 

Up to this point we have considered only coarse delay estimation by fractional 
testdelays, which is clearly enough to achieve an estimation error less that half 
a chip duration. The accurary of the estimate can then be achieved in tracking 
mode. Of course, accuracy could be achieved also by using smaller stepsizes for the 
testdelays. However, differential correlations enable actually solving the fractional 
part of the delay. Namely, we show in the following that in a frequency-flat fading 
channel it is enough to solve a system of two second-order polynomials for each chip 
interval. 

Assume that only one resolvable path exists. Suppose d+6is the propagation delay, 
so that the differential correlation matrix becomes R(r) = a(r)cucfj. Recall that 
cii = Cn + cii by Eq. (20). According to Eqs. (4) and (5), we have moreover 
cn = (1 - S)gi{d) + 5gi{d -f 1), where gi(d) and gi(d + 1) are defined in Eq. 
(24), being just a d and d + I times shifted rephcas of the desired chip sequence, 
respectively. Suppose now that we use testcodes exatly gi(d) and gi(d+ 1). This is 
because the true delay lies between d and d + l. Considering the MF-type approach, 
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we immediately have (see details in the Appendix B) 

|gi(dfR(T)gi(d)| = \a{r)\[{l - S)m + Svil)? (30) 
|gi(d+lfR(T)gi(d + l)| = |a(r)|((l - + ^77(0)l^ (31) 

where rjii) gi(rf)^gi(rf + »)• Notice that these quantities are known, since they 
are (up to scaling) just values of the code autocorrelation function with time lag i 
chips. For example, for m-sequences 77(0) = C, and t/(1) = -1. Hence, only two 
unknowns, |d(r)| and 6 are included in Eqs. (30) and (31). The solutions for 5 can 
be found (see details in the Appendix C) to be equal to 

*^^^^-(r;(0)-;7(i))(i±W)' ^ ^ 

where 

p{d) = 



gx{d)TR{r)g,{d)\ (33) 



gi(d+l)m(r)gi(d+l)r 
Eq. (32), which we label ADC-MF (Accurate DC-MF), is a closed form expression 
for the timing offset S, for which we need only the information about the hypotheses 
d and d 4- 1. ADC-MF is summarized in Table 1. 

Also MUSIC-type approach could be used. This is because R(r) = a(r)cuc?'i has 
only one non-zero eigenvalue, and the related eigenvector is Cu/||cu||. Replacing 
R(r) by U(t)U(t)^ = Ciicfi/||cn|P is Eqs. (30) and (31), we would have 

|gi((i)^U(r)U(r)^gi(d)| = jj-L_[(l - <J)7,(0) + <Jr?(l)l2 (34) 
|g,(d+l)^U(r)U(r)^gi(d+l)l = jj^Kl - %(1) + MO)f . (35) 

Since ||cii|| is unknown, as being a function of <J, we end up with Eqs. (32) and (33). 
Naturarly, R(t) is replaced by U(r)U(r)^ in Eq. (33). Analogously, we label this 
approach as ADC-MUSIC (Accurate DC-MUSIC), and summarize it also in Table 
1. 
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Remark 4. A slightly modyfied version of Eq. (32) could be used also in case of 

multipaths. This is because now 

— L 
R-C^) = X^a<(')cucfi , i.e. Eq. (19) 

= i:a,(r)[(l-M'giW)giWf 

1=1 

+ 25,(1 - Si)gi{di)gM + If + ShM + l)giidi + If]. (36) 

Letting now d = rfi, equations (30) and (31) would accordingly have more general 
forms 

gi{d,fR{T)gM = a,(T)lil ^ 6^U0) + St,{1)]' 

gi(di + l)^R(r)gi(di + l) = ai(r){(l-(J)77(l) + (57,(0)p 

+ E ^'(t)[(1 - Si)r}{di ~ d/ - 1) + (Jr;(rfi - di)p8) 

1=2 

Assuming impulse-like code autocorrelation, i.e. T){i) « 0 for |i| > 1, we would end 
up with (30) and (31), where in addition 7/(1) = 0. This would accordingly lead to 
an estimator 

'=TTm- '''' 

3. NUMERICAL EXPERIMENTS 

We compared the methods in the downlink environment with Rayleigh fading chan- 
nel. Data rate was assumed to be fd = 16 kbit/s, carrier frequency fc = 1.8 GHz, 
and mobile speed i; = 50 km/h. This results in normalized Doppler shifts at most 
~^fc/fd = 5.21 • 10"*^. Typical fading process at the symbol level is seen in Figure 

''light 

3. 



The delays were estimated in parallel manner, and the strongest peak from the delay 
sprectum was selected. Acquisition was declared if this peak corresponded to the 
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strongest path. Hence, our requirement is more stringent than the other common 
definition, for which it is sufficient just to find any multipath component. The 
measured"performance parameters were probability of acquisition, mean acquisition 
time, and root-mean-square-error (RMSE). Mean acquisition time (in symbols) was 
measured as [23) 



T^CQ = "^""T ^"^^^^^ Te, (40) 



where Tp and Tc are the penalty time, in case of misacquisition and chip duration, 
respectively. In our simulations we had Tp = 500 symbols and Tc = ^ = ^ symbols. 
PrACQ is the probability of acquisition. The other performance parameter, RMSE, 
measures the accuracy of the delay estimate according to 



RMSE = 

where 5 is the number of simulations. 



1 ' 



M5E(d-d)^ (41) 



The methods for delay estimation were: DC-MF (23), DC-MUSIC (28), ADC-MF 
and ADC-MUSIC (Table 1), MUSIC (29), CMOE [17] , and conventional non- 
coherent and coherent MF [4]. Any parameters such as time- varying path strengths, 
and codes of the interfering users were not assumed to be known. Only the code of 
the desired user was known. In addition, a preamble = 1 was available for the 
desired user A; = 1. Gold codes of length C = 31 were used. Signal- to-Noise Ratio 
(SNR) in the chip-matched filter output was always 10 dB. In frequency-selective 
fading channel the number of paths was L = 3 if not otherwise stated. Relative 
differences in paths powers were 5 dB, and SNR was always set with respect to 
the strongest path. The number of users K was varied from if = 5 to A' = 25. 
Therefore, the total number of sources was 2KL = 6K = 30, . , . 150, while the 
maximum signal subspace dimension was C = 31. Therefore, highly loaded systems 
arise. We noticed that the performance of MUSIC became better by selecting the 
estimated signal subspace dimension 2K instead of 2KL. This can be justified also 
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by the fact that usually it is impossible to know the number of paths, while number 
of users can be known by uplink-downlink communication. In one simulation, each 
parameter value was fixed, and the delays were estimated from the peaks of the 
delay spectra. All the measured quantitites are averaged over 1000 simulations. 

In the first experiment we studied the eflfect of increasing the length of the preamble. 
The system included K = 20 users and the desired user had either equal or 5 dB 
weaker signal power than the others. From Figure 4 we see that non-coherent MF 
failed in near-far scenario, which was predictable, as well as MUSIC due to high 
system load. The zero mean fading process also seemed to be too fast for the co- 
herent MF, even though it was in this setting the most robust against MAI. CMOE 
performed better the more samples were available. This was because the estimate 
for the inverse of the autocorrelation matrix R(0) thus became more accurate. How- 
ever, in the presence of MAI, see Figure 5, CMOE did not really managed to catch 
at least the strongest path first. The fact that the proposed methods are asymp- 
totically invariant with respect to the interfering users and noise can be predicted 
from Figures 4 and 5, in which the probability of acquisition seems to tend to the 
to 1, whether the MAI is 0 or 5. In this experiment, M = 500 symbols were enough 
for them to achieve almost the same performance whether the interfering users were 
5 dB stronger or not. The figures also tell the fact that the performance of the 
proposed methods depends essentially on the lenght of the preamble, because the 
goal is to average interference out. For this reason, it is clear that the methods fail 
if the preamble is too short. On the other hand, although coherent MF was the best 
among the others when short preamble is used, its performance remained at mod- 
erate level. In Figures 6 and 7, the mean acquisition times in symbols are plotted. 
Due to the reliable estimation, the acquisition time for the proposed methods was 
usually no more than the length of the preamble. As can be seen from Figure 5, 
coherent MF was in average the fastest if only a short preamble was available. One 
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can also notice, that DC-MF seemed to have better acquisition cabability of the 
strongest path than DC-MUSIC. The reason for this is that MF-type approaches 
work ingeneral well as the code autocorrelation function is impulse-like. To show the 
power of DC-MUSIC when the code autocorrelation is not that spiky, we assigned 
K = 2 users random codes of lenght C = 7 in L = 3 path Rayleigh fading channel. 
The paths had 5 dB differences in power. Figure 8 shows the average number of ac- 
quired paths as a function of preamble length. It is seen that already with preamble 
lenght M = 100 DC-MUSIC found the most of the paths. Table 2 shows the case 
M - 500 in more detail. Namely, it tells the number of cases when each path was 
acquired. It is quite noticeable how frequently DC-MUSIC acquired also the weaker 
paths (-5 dB and -10 dB paths with probabilities 0.917 and 0.629, respectively) 
compared to the other methods, which had serious problems in acquiring even the 
second strongest path. 

In the second experiment we studied the effect of system loading, i.e. when the 
number of users was increased. The system included K = 5to K = 25 users with 
MAI equal to 0 dB. The length of the preamble was M = 500 symbols. Figure 9 
shows the achieved probability of acquisition. Although the true model order (equal 
to 2KL) was quite high, the model order underestimation (estimated as 2K) helped 
MUSIC in the case of A" = 5 users. This was not the case for the other system 
loads. Coherent MF had nearly the same performance regardless of the system 
load, althoug the performance level again was moderate. The difference of the other 
remaining methods (DC-MF,DC-MUSIC,CMOE, and MF) in this experiment was 
too small up to i<' = 15 users to draw conclusions. This was mainly because of the 
MAI free case. On the other hand, especially in the case K = 25, which means quite 
a heavy system load, the proposed methods clearly performed the best. It is also 
seen that coherent MF, together with the proposed methods are the most robust 
against system loading. But among those only the proposed ones can exploit the 
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time-correlatedness of the fading process, resulting with better performance level. 
Figure 10 tells the corresponding mean acquisition times. Assuming the bit rate of 
16 kbit/s; these times for the case = 25 would equal the following: DC-MF and 
DC-MUSIC 33 ms, coherent MF 44 ms, CMOE 52 ms, and non-coherent MF 54 
ms. Thus, at least 25% improvement would be gained by the proposed methods, 
compared to any reference ones. 

In the third experiment we studied the behavior of the methods when increasing 
the number of inter-symboHnterference (ISI). The ISI was generated by increasing 
the number of paths from L = 2 to L = 10. One of the paths was always 5 dB 
stronger, which was intended to be acquired. System included K = 20 users, and 
M = 500 symbols were used again as a preamble. MUSIC was no more applied due 
to its obvious failure in highly loaded systems (see the first experiment). Hence, 
the only method which would need to know or estimate the number of paths L, is 
DC-MUSIC. Recall that DC-MUSIC computes L eigenvectors. In this experiment 
DC-MUSIC always assumed. L = 3. Figure 11 shows that proposed methods clearly 
resist ISI the best. Moreover, DC-MUSIC seems not to be that sensitive to the 
model order mismatch. 

In the fourth experiment, we studied the achieved accuracy in a system with K = 20 
users. Although MUSIC is a high-resolution method, it was not applied, again, due 
to its failure in this setting. The channel was frequency-flat fading, and ADC- 
MF and ADC-MUSIC were applied for the first time. The accuracy for all the 
other methods was achieved by using more delay candidates. More precisely, the 
increment for the delay candidate was always as small as chip. Recall that ADC- 
MF and ADC-MUSIC uses an increment equal to one chip, only. Figure 12 shows 
the achieved accuracies in RMSE. It can be stated that the high-resolution feature 
of MUSIC is now seen in the performance of ADC-MUSIC, which reached an RMSE 
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of only 0.02 - 0.04 chips, depending on the system load. CMOE seems also to yield 
quite accurate. estimates, and performs better than ADC-MF, which reached RMSE 
of 0.04 — U.08 chips. However, it must be pointed out that there was no MAI in this 
experiment. It is hence expected that in the presence of MAI, the RMSE of CMOE 
as well as coherent and non-coherent MF will degrade significantly, as did in the 
first experiment. 

The final experiment was dedicated for the comparison of the accuracy of ADC-MF 
and ADC-MUSIC to the Cramer-Rao lower bound (CRB). Although the estimators 
turned out to be approximately unbiased, the bias was estimated for each fixed 
parameter value. This was taken into account in the final RMSE calculation, giving 
an adequate measure for CRB comparison. The system included jFf = 5 users, and 
the average SNR was 10 dB. Single-path time-invariant channel was assumed, and 
hence the data model coincide e.g. that of [10], where the CRB was derived [pp. 
88-89], The RMSE curves are seen in Figure 13, from which we notice that neither 
methods attain the CRB. This was the case also in [9] with MUSIC. However, it 
must be pointed out that the preamble is not yet optimally utilized in ADC-MF 
nor in ADC-MUSIC. This is because the contribution of interference and noise was 
tried to suppress by estimating only one differential correlation matrix R(r), with 
r = 2 in all the experiments* But especially with slowly fading channels, strong 
time correlatations takes place with much larger time intervals, which suggests the 
possibility of achieving even more interference suppression, and consequently better 
accuracy. 



DERI\^ATION OF EQUATION (19) 
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First, from Eq. (8) we have 



GA(t) = (£ii,Cn,...,£ii,Cii„...,£|fir,C/fL] 



/i?ai(T)l 



(42) 



= (/i?Q!i{T)(cii + cu) • • -^iaiMCcii + cil) 0- . -0] 
= [/ufai(r)cii---/iiai(r)cii 0---0], 



(43) 
(44) 



1 1 
1 1 



where the second equation follows from 1 = 
definition = fi^ 4-Cij. Therefore we have 

L 

GA(r)G^ = Y,^^fal{T){cu&l + Ci,c[,), 



, and the third is due to 



(45) 



from which (19) immediately follows. 

DERIVATION OF EQUATIONS (30) AND (31) 

Recall from section 3.3, that 

R(r) = Q:(r)cnC^i , and 
cn = (l-<J)gi(rf) + 5gi(d + l). 

Thus, we first may write 

cucl = (l-<5)'gi(rf)gi(rf)^ + 5'gi(rf+l)gi(d+l)^ 

+ 5{1 - <J)gi(d)gi(d + 1)^ + <y(l - S)gi{d + l)gi(rf)^. 

Defining 7j(i) = gi(d)''gi(d + i), we have 



(46) 

(47) 



(48) 



gx(d)^R(r)gi(d) = a(r)[l - 6)M0? + 6^1? + S{1 - 5)^(0)77(1) + 6{1 - SUlUQ)] 
= a{T)[{l-6)viO)+6ri{l)?, (49) 
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from which Eq. (30) follows. The derivation of Eq. (31) is entirely similar and is 
omitted. 

DERIVATION OF EQUATION (32) 

To shorten notations, define A = gi(d)^R(r)gi(ci) and B = gi(d+l)^R(r)gi(d+l). 
Taking squareroots on both sides of Eqs. (30) and (31), we have 



which is a system of equations with two unknowns, 5, and yja{r). Dividing Eq. 
(50) by Eq. (51) we obtain 




(50) 




(51) 





(52) 



By simple algebraic manipulations it leads to 




(53) 



from which Eq. (32) is obtained by denoting P{d) = 
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According to a further aspect of the invention, a synchronization method in a spread 
spectrum receiver is provided, in which method matched filtering is employed and 
the delay of received signal is estimated. In an advantageous embodiment of the 
invention, in matched filtering, a set of test delays is scanned through with the 
5 pseudonoise code used in matched filtering, consecutive test delays having a delay 
difference of one chip interval, and each different test delay producing a sequence 
of matched filtering results, and 

each said sequence of matched filtering results is correlated with a delayed copy of 
each said sequence of matched filtering results for producing a correlation result 
10 value corresponding to the test delay producing said sequence of matched filtering 
results, 

and in that, that delay estimation comprises steps, in which 

- two further test delays are calculated from each pair of two consecutive correlation 
result values, 

15 - it is checked for each of said two further test delays, if the test delay is between the 
test delay values corresponding to said pair of two consecutive correlation result 
values, and if it is, the test delay is added to the set of test delays, and 

- a delay from the set of test delays is selected to be a result of delay estimation. 

20 In the step of calculation of two further test delays, preferably equation (32) or 
equation (39) is used. 

The present aspect of the invention is illustrated further in figure 14, which shows a 
part of a block diagram of a spread spectrum receiver. Figure 14 shows a matched 

25 filter block 10, a delay block 15, a correlation block 20, and a delay estimation block 
30. In the figure, r(t) denotes the stream of sample vectors, and PN denotes the 
pseudonoise i.e. code fed into the matched filter block 10. Different matched 
filtering results are obtained using a set of test delay values for the pseudonoise code 
PN fed into the matched filtering block, and each stream of matched filtering results 

30 is correlated with a delayed copy of itself for producing a correlation result value. In 
the delay estimation block, further test delays are calculated on the basis of said 
correlation result values, and the test delay producing the best niatch to the received 
samples is selected as the estimated delay. It is also possible to select a plurality of 
the test delays, if the signal from more than one paths is desired to be detected. 

35 

The selection of a test delay to be a result of the delay estimation can 
advantageously be performed using condition (26). 
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In one advantageous embodiment of the invention, more than one delays are 
estimated. In such an embodiment, after the first test delay is selected as a result of 
the estimation, the test delay is removed from the set of test delays as well as all 
those test delays, whose delay difference to the selected and removed test delay is 
5 less than one chip interval. The next delay estimate is then selected from the 
remaining set. 

According to a further aspect of the invention, a synchronization method in a spread 
spectrum receiver is provided, in which method signal subspace estimation is 

10 employed. According to an advantageous embodiment of the invention, sampled 
data and at least one delayed version of sampled data is used in the signal subspace 
estimation in the method, and the method further comprises the step of calculation 
of a signal subspace which is substantially the same as a signal subspace calculated 
of a combination of at least two matrices, which matrices have the form of MM^, 

15 where M is a matrix composed of a first set of data samples, and is the 
hermitian transpose of a matrix composed of a second set of data samples, said 
second set of data samples having a time difference of n samphng intervals to 
samples of said first set of data samples. 

20 The time difference n can be a positive or a negative integer, whereby the second set 
of data samples may be an earlier set or a later set as the first set of data samples. 
We note here that with small values of n, i.e. when the absolute value of n is smaller 
than the size of the sets of data samples forming the matrices M and M„, the 
matrices may have some overlap. In other words, the matrices M and M„ are 

25 formed of sample vectors r^...r„ and r„+„...r^+„, respectively, where is the first 
sample vector in M and the last sample vector in M . 

In various embodiments of the invention, the signal subspace can be calculated in 
many different ways, for example directly from the data vectors. The signal 
30 subspace can also be calculated from the combination of MM^ matrix products. In 

a further advantageous embodiment of the invention, in the step of calculation of a 
signal subspace, said signal subspace is calculated of a combination of at least two 
matrices, which matrices have the form of MM}jj , where M is a matrix composed 
of a first set of data samples, and mH is the hermitian transpose of a matrix 
35 composed of a second set of data samples, said second set of data samples having a 
time difference of n sampling intervals to samples of said first set of data samples. 
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Said combination is preferably a linear combination, which is useful for 
symmetrization purposes. However, the invention is in no way limited only to a 
linear combination of matrices of the form of MMjf . 

5 Calculation of a signal subspace of matrix A refers to the calculation of such a 
subspace, whose basis vectors span the same subspace as the eigenvectors of matrix 
A corresponding to the largest eigenvalues of matrix A. Calculation of signal 
subspace is in itself well known to a man skilled in the art, wherefore it is not 
described in more detail in this application. 

10 

The present aspect of the invention is illustrated further in figure 15. The figure 
shows a part of the block diagram of a spread spectrum receiver. Figure 15 shows a 
signal subspace estimation block 60, whose function is to produce an estimate of the 
signal subspace for use in delay estimation. According to the invention, delayed 
15 sample matrices are used the signal subspace estimation in addition to sample 
matrices collected without delay. This is represented by delay blocks 50 in figure 15. 
As the figure illustrates, there may be more than one delay, which is to say that the 
signal subspace estimation can be performed on the basis of a plurality of products 
MM^ , where each delay corresponds to a different value of n. As stated before, the 

20 length of the delay is an integer multiple of the sampling interval. In the figure, r(t) 
denotes the stream of sample vectors. 

According to an advantageous embodiment of the invention, a synchronization 
methods according to an advantageous embodiment of the invention is used in a 
25 mobile conraiunication means of a cellular teleconmiunication system. 

In an advantageous embodiment of the invention, a mobile conmiunication means is 
provided, which mobile communication means comprises 

- a subspace estimation block arranged to use sampled data and at least one delayed 
30 version of sampled data is used in the signal subspace estimation, and 

- means for calculation of a signal subspace which is substantially the same as a 
signal subspace calculated of a combination of at least two matrices, 

which matrices have the form of MM^ , 

where M is a matrix composed of a first set of data samples, and is the 
35 hermitian transpose of a matrix composed of a second set of data samples, said 
second set of data samples having a time difference of n sampling intervals to 
samples of said first set of data samples. 
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According to an advantageous embodiment of the invention, a synchronization 
methods according to an advantageous embodiment of the invention is used in used 
in a global positioning system (GPS) receiver. 

5 In an advantageous embodiment of the invention, global positioning system receiver 
is provided, which receiver comprises 

- a subspace estimation block arranged to use sampled data and at least one delayed 
version of sampled data is used in the signal subspace estimation, and 

- means for calculation of a signal subspace which is substantially the same as a 
10 signal subspace calculated of a combination of at least two matrices, 

which matrices have the form of MMj' , 

where M is a matrix composed of a first set of data samples, and Mj; is the 
hermitian transpose of a matrix composed of a second set of data samples, said 
second set of data samples having a time difference of n sampling intervals to 
1 5 samples of said first set of data samples. 

According to an advantageous embodiment of the invention, a synchronization 
methods according to an advantageous embodiment of the invention is used in a 
receiver of a wireless local area network (WLAN) terminal. 

20 

In an advantageous embodiment of the invention, a terminal of a wireless local area 
network is provided, which terminal comprises 

- a subspace estimation block arranged to use sampled data and at least one delayed 
version of sampled data is used in the signal subspace estimation, and 

25 - means for calculation of a signal subspace which is substantially the same as a 
signal subspace calculated of a combination of at least two matrices, 
which matrices have the form of MMU , 

where M is a matrix composed of a first set of data samples, and is the 
hermitian transpose of a matrix composed of a second set of data samples, said 
30 second set of data samples having a time difference of n sampling intervals to 
samples of said first set of data samples. 

According to an advantageous embodiment of the invention, the synchronization 
methods according to various embodiments of the invention are used in radio link 
35 receivers employing spread spectrum technology. For example, any synchronization 
method described in this specification can be used for example in the well-known 
Bluetooth device. 
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The present invention provides an improved method for performing initial 
synchronization. The invention can be used with any tracking mechanisms and other 
mechanisms associated with spread spectrum reception, wherefore the invention is 
not hmited to to being used with any specific tracking mechanisms or other 
5 mechanisms associated with spread spectrum reception. 

The equations described in this specification are meant to be examples only, since 
the same underlying ideas can be expressed in a multitude of different ways using 
different mathematical phraseology, as a man skilled in the art very well knows and 
10 understands. Therefore the invention is not limited to exactly those equations 
described in this specification nor to any other specific set of equations. 

The inventive method is applicable in spread spectrum receivers. Practical examples 
of such applications which are commercially important at the time of writing this 
15 patent application are mobile stations of WCDMA cellular systems such as the 
UMTS, IS-95, and IMT-2000 cellular systems, global positioning system (GPS) 
receivers and wireless local area network (WLAN) terminals. 

In view of the foregoing description it will be evident to a person skilled in the art 
20 that various modifications may be made within the scope of the invention. While a 
preferred embodiment of the invention has been described in detail, it should be 
apparent that many modifications and variations thereto are possible, all of which 
fall within the true spirit and scope of the invention. 
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Denote Q the set of delay candidates. 

1. Estimate the differential correlation matrix with time lag r from available 
vector samples as R(r) = ^ Em^i ^mr^^r- 

2. Let U(t) contain L principal eigenvectors of R(t) + R(r)^, where L is the 
number of paths to be acquired. 

3. Set fi = {0, 1, . . . . C}, i.e. multiples of a chip duration. 

4. For each deil- {C} do: 

(a) Compute two candidates for fractional part of the delay, according 
to Eqs. (32) and (33). In the latter equation, replace 

• R(r) by R(r) in case of ADC-MF, 

• R(r) by U(r)U(r)^ in case of ADC-MUSIC 

(b) Insert d + (J(d) to if 6{d) 6 (0, 1). 

5. Obtain one delay estimate by choosing the strongest peak from the delay 
spectrum 

|cx(x)^ec,(x)| I e = R(r), ADCMF ^^^^ 
\MxW I e = U(r)U(r)^, ADC-MUSIC 

where Ci(x) is a testcode of form Eq. (24), corresponding to delay x ^ 

6. If other paths are acquired, too, choose such a next strongest peak of Eq. (54), 
which does not correspond to a delay closer than one chip to any already found 
delays. 



Table 1: Algorithms ADC-MF and ADC-MUSIC for timing offset estimation. 
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path 


1« 




3rd 




power 


0 dB 


—0 dB 


-10 dB 




delay (chips) 


1 + di 










#/1000 


#/1000 


#/1000 


E/3000 


DC-MUSIC 


969 


917 


629 


2515 


DC-MF 


975 


710 


414 


2099 


MF 


958 


620 


357 


1935 


coh. MF 


878 


570 


404 


1851 


CMOE 


953 


765 


433 


2151 



Table 2: Number of cases a certain path was acquired in a -ff = 2 user system with 
random codes of length C = 7. Preamble length was M = 500 symbols, and average 
SNR was 10 dB. The fractional parts <Jf, i = 1,2,3 were uniformly distributed over 
(0,1). 
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Claims 



1. Synchronization method in a spread spectrum receiver in which method signal 
subspace estimation is employed, characterized in that 
5 sampled data and at least one delayed version of sampled data is used in the signal 
subspace estimation, and in that, that the method further comprises the step of 
calculation of a signal subspace which is substantially the same as a signal subspace 
calculated of a combination of at least two matrices, 
which matrices have the form of MM^ , 
10 where M is a matrix composed of a first set of data samples, and is the 
hermitian transpose of a matrix composed of a second set of data samples, said 
second set of data samples having a time difference of n sampling intervals to 
samples of said first set of data samples. 

15 2. A synchronization method according to claim 1, characterized in that 

in the step of calculation of a signal subspace, said signal subspace is calculated 
of a combination of at least two matrices, 
which matrices have the form of MM" , 

where M is a matrix composed of a first set of data samples, and is the 
20 hermitian transpose of a matrix composed of a second set of data samples, said 
. second set of data samples having a time difference of n sampling intervals tO; 
samples of said first set of data samples. 

3. Synchronization method in a spread spectrum receiver in which method matched 
25 filtering is employed and the delay of received signal is estimated, characterized in 
that in the method, 

in matched filtering, a set of test delays is scanned through with the pseudonoise 
code used in matched filtering, consecutive test delays having a delay difference of 
one chip interval, and each different test delay producing a sequence of matched 
30 filtering results, and 

each said sequence of matched filtering results is correlated with a delayed copy of 
each said sequence of matched filtering results for producing a correlation result 
value corresponding to the test delay producing said sequence of matched filtering 
results, 

35 and in that, that delay estimation comprises steps, in which 

- two further test delays are calculated from each pair of two consecutive correlation 
result values, 
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. it is checked for each of said two further test delays, if the test delay is between the 
test delay values corresponding to said pair of two consecutive correlation result 
values, and if it is, the test delay is added to the set of test delays, and 

- a delay from the set of test delays is selected to be a result of delay estimation. 

5 

4. A mobile station of a cellular telecommunication system utilizing spread 
spectrum communication, 

characterized in that it comprises 

- a subspace estimation block arranged to use sampled data and at least one delayed 
10 version of sampled data is used in the signal subspace estimation, and 

- means for calculation of a signal subspace which is substantially the same as a 
signal subspace calculated of a combination of at least two matrices, 

which matrices have the form of M . 

where M is a matrix composed of a first set of data samples, and Mj^ is the 
15 hermitian transpose of a matrix composed of a second set of data samples, said 
second set of data samples having a time difference of n sampling intervals to 
samples of said first set of data samples. 

5. A global positioning system receiver, 
20 characterized in that it comprises 

- a subspace estimation block arranged to use sampled data and at least one delayed 
version of sampled data is used in the signal subspace estimation, and 

- means for calculation of a signal subspace which is substantially the same as a 
signal subspace calculated of a combination of at least two matrices, 

25 which matrices have the form of MM^ , 

where M is a matrix composed of a first set of data samples, and is the 
hennitian transpose of a matrix composed of a second set of data samples, said 
second set of data samples having a time difference of n samphng intervals to 
samples of said first set of data samples. 

30 
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6. A tenninal of a wireless data communications network, 
characterized in that it comprises 

- a subspace estimation block arranged to use sampled data and at least one delayed 
version of sampled data is used in the signal subspace estimation, and 
5 - means for calculation of a signal subspace which is substantially the same as a 
signal subspace calculated of a combination of at least two matrices, 
which matrices have the forai of MM" , 

where M is a matrix composed of a first set of data samples, and M" is the 
herraitian transpose of a matrix coniposed of a second set of data samples, said 
10 second set of data samples having a time difference of n sampling intervals to 
samples of said first set of data samples. 
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